Infrared and ultraviolet properties of the Landau gauge quark propagator 
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We present a current summary of a program to study the quark propagator using lattice QCD. We use the 
Overlap and "Asqtad" quark actions on a number of lattice ensembles to assess systematic errors. We comment on 
the place of this work amongst studies of QCD Green's functions in other formulations. A preliminary calculation 
of the running quark mass is presented. 



1. INTRODUCTION 

The quark propagator is a building block of 
Quantum Chromodynamics and has been stud- 
ied in numerous formulations |1I2I3| . The lattice 
regularisation provides a method for a full non- 
perturbative calculation from first principles. In 
the infrared, the quark propagator displays dy- 
namical mass generation, but in the ultraviolet it 
takes its perturbative, running form. 

By computing the quark propagator on the 
lattice we can connect with perturbation the- 
ory and with many models. We can also con- 
tribute to our understanding of the inner work- 
ings of the lattice simulations themselves. It has 
been seen, however, that detailed lattice studies 
of the quark propagator require sophisticated ac- 
tions |4l5lfil7l8l§| . 

We present some results for the quark propa- 
gator in Landau gauge from a variety of lattice 
spacings, two different physical volumes and two 
different quark actions. Specifically we use an im- 
proved Staggered action, "Asqtad," ^U] and the 
Overlap quark action. Results from Asqtad and 
Overlap actions are in good agreement. We will 
compare our results to some model calculations 
and compare our chiral extrapolation with recent 
Dyson-Schwinger equation predictions. Not only 
do we see stable infrared results, but the ultra- 
violet data is good enough for us to attempt es- 
timates of the chiral condensate and the running 
quark mass. A much more detailed presentation 
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can be found in Ref. jllj . The quark propagator 
was also studied in Laplacian gauge in Ref. UJ. 

2. THE LATTICE QUARK PROPAGA- 
TOR 

In the (Euclidean) continuum, Lorentz invari- 
ance allows us to decompose the full quark prop- 
agator into Dirac vector and scalar pieces 



S{p 2 ;m) 
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Asymptotic freedom means that at large momen- 
tum transfer, p 2 — > oo, the quark propagator ap- 
proaches its free form, 
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where m is a current quark mass. On the lattice 
in the free case (U^x) — 1), the quark propagator 
takes the form, 
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where the function q^ipp.) depends on the quark 
action. For the Asqtad action it is 0, 



= sin(p M ) [1 + - sin 2 (p Al )] , 
and for the Overlap action it is 8 , 



q M = 2m, 
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Figure 1. Tree-level behaviour: q vs. p for the 
Overlap (top curve) and Asqtad (bottom curve) 
quark actions. This is the inverse of the propaga- 
tor for a massless quark. The diagonal line is the 
continuum case. 



where 

A = -m w + — 

% = sin(p M ) 

kf, = 2sin(p Al /2), 

and m w is the regulator mass. These functions 
are plotted in Fig. ^ 

The difference between Eq. and Eq. © is a 
pure lattice artefact. For the Dirac vector piece 
this can be trivially corrected for by considering 
it to be a function not of p^, but of q^. This 
only affects the ultraviolet behaviour of Z, and 
asymptotic freedom means that we may hope that 
this correction will persist in the interacting case. 
This choice of "kinematical momentum" has been 
seen to be efficacious in the removal of hypercu- 
bic artefacts jo] and the philosophy is the same 
as that pursued with great success in studies of 
the gluon propagator ^21- Tree- level correction 
does not prescribe such a procedure for the Dirac 
scalar piece, M, so we shall adopt an empirical 
approach. 



Table 1 



Simulation parameters 



Action 


a 


a a (fm) 


L a 


x T 


Asqtad 


4.38 


0.167 


16 3 


x 32 


Asqtad 


4.60 


0.124 


12 3 


x 24 


Asqtad 


4.60 


0.124 


16 3 


x 32 


Overlap 


4.286 


0.190 


8 3 


x 16 


Overlap 


4.60 


0.124 


12 3 


x 24 


Overlap 


4.80 


0.093 


16 3 


x 32 



3. RESULTS 

3.1. Lattice artefacts 

The simulation parameters are given in Tabled 
In all cases an 0{a 2 ) Symanzik improved gluon 
action was used. The propagators were rotated to 
0(a 2 )-improved Landau gauge JT3J by standard 
methods. 

As we have both Overlap and Asqtad results 
for one set of lattices we can make a direct com- 
parison. This is shown in Fig. [2] for one choice of 
comparable quark masses. Z is shown as a func- 
tion of q - chosen appropriately for each action 
- and there is complete agreement, within errors. 
The Asqtad action has notably better rotational 
symmetry, presumably reflecting the removal of 
0(a 2 ) errors. Recently, one group has attempted 
to correct the poor rotational symmetry of the 
Overlap action by extrapolating in the hypercu- 
bic artefacts The quark mass function for 

the Asqtad action is shown as a function of q, 
which has been seen to produce slightly better 
results jjj , but we use the standard lattice mo- 
mentum, p for the Overlap action. This choice 
produces the best agreement between the two ac- 
tions. 

Figure [3] shows the Asqtad quark propagator 
on two different size lattices with the same lattice 
spacing. The smaller lattice has a spatial length 
of about 1.5 fm, the larger about 2.0 fm. The 
mass functions lie on top of each other, indicat- 
ing that finite volume effects are small, if present. 
For Z however, there is a clear difference in the 
low momentum region. Infrared suppression of Z 
is a long-standing prediction of Dyson-Schwinger 
equation studies and the larger lattice results are 
in good agreement with recent DSE results |15|. 
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Figure 2. Comparison of the actions: Quark mass 
function (top) and Z function (bottom) for mass 
m ~ 80 MeV at j3 = 4.60. These calculations 
were performed on a 12 3 x 24 lattice. 



Figure 3. Finite volume effects: Asqtad quark 
mass function (top) and Z function (bottom) for 
mass m ~ 57 MeV at p = 4.60. Comparison on 
12 3 x 24 lattice (open circles) and 16 x 32 lattice 
(solid squares). 



but at least some of the suppression may be due 
to the finite volume. 

As we have the Overlap quark propagator at 
three lattice spacings we can examine its scaling 
behaviour. This is shown in Fig. For both func- 
tions we see good consistency across these lattice 
spacings. Once again, we have plotted Z against 
q and M against p. This gives the best agreement 
between results at different lattice spacings. The 
violation of rotational symmetry seen in Z - quite 
strong on the coarsest lattice - becomes smaller 
with the lattice spacing, but the overall form is 
stable. The infrared behaviour of the mass func- 
tion is almost unaffected by the change in lattice 
spacing, but the ultraviolet shows some (not un- 
expected) sensitivity. 



The quark propagator is, of course, a function 
of the bare quark mass, and this dependence is 
illustrated in Fig. For the Z function we show 
results for three quark masses from the Asqtad 
action, and indeed some mass dependence is ob- 
served in the infrared. The heaviest quark has 
about four times the mass of the lightest one, so 
this dependence is rather weak. We expect the 
ultraviolet to be insensitive to the quark mass 
and this is in fact observed. The ultraviolet tail 
of the mass function is proportional to the bare 
quark mass and this is observed in Fig. using 
the Overlap action. As the mass increases the in- 
frared enhancement becomes less significant. It 
has already been remarked elsewhere [7] that for 
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Figure 4. Scaling: Quark mass (top) and Z (bot- 
tom) for the Overlap action at three lattice spac- 
ings. Bare quark mass m ~ 106 MeV. Z is shown 
as a function of q while M is a function of p. This 
is seen to provide the best agreement between re- 
sults on different lattice spacings. 



Figure 5. Mass dependence: Asqtad Z function at 
j3 = 4.60, 16 3 x 32, for bare quark masses m ~ 57, 
115 and 229 MeV. There is some weak mass de- 
pendence in the infrared. Overlap mass function 
at (3 = 4.60, 12 3 x 24, for bare quark masses 
to ~ 106 - 531 MeV. For large momentum, M(p 2 ) 
is proportional to the quark mass and as the mass 
increases, M(p 2 ) becomes flatter. 



small quark masses the infrared part of the mass 
function is insensitive to the quark mass, but 
at around the strange quark mass the shape of 
the mass starts to change. In the non-relativistic 
limit M(p 2 ) would be a constant. 

3.2. Comparison with other methods 

We extrapolate the quark mass function to the 
chiral limit using a quadratic fit. One reason for 
this is that it is a convenient quantity to compare 
with other calculations. In Ref. the authors 
suggested a novel form of the quark mass function 
which actually vanishes at zero four-momentum. 



The ansatz 



M(q 2 ) = 



where 



Q(q 2 



l-y/l-g 2 Q 2 (g 2 
Q(<f) 



-(ilnp) 16?r 2 
~2N C F~ 



-9 2 /A 2 



(6) 



(7) 



based on the Constrained Instanton Model and 
QCD sum rules with non-local condensates, is 
one of the curves shown in Fig. This form 
is not favoured by the lattice data. It is fas- 
cinating that such a form should be consistent 
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Figure 6. Comparison with two models: The 
mass function that vanishes at zero four- 
momentum is from a model of Dorokhov and Bro- 
niowski [3], the other is from an instanton calcu- 
lation by Diakonov [2]. 



with the phenomenological constraints discussed 
in Ref. |3] and demonstrates the utility of first- 
principles calculations. 

Figure also shows a regular instanton model 
calculation from Ref. There the quark mass 
function takes the form 



M(q) = M(0)F 2 (z) 
F(z) = 2z(l (z)K 1 (z) - h(z)K (z) 



(8) 



h{z) 



P 

Z = q 2 



where the I's and K's are the modified Bessel 
functions. The parameters M(0) — 350 MeV and 
p^ 1 = 600 MeV are derived from the standard 
values of the instanton ensemble. This two pa- 
rameter form has a good correspondance with the 
lattice data. 

The interaction between lattice and DSE stud- 
ies has been particularly active recently. See 
Refs. |15I16I17I18| for a few examples. Amongst 
other things, this is providing valuable insight 
about the validity of various DSE truncation 
schemes. In return, DSE calculations can be done 
at quark masses and momenta difficult to access 
on the lattice and can also provide information 



on the time-like region. 

In Ref. \R>\ there is a comparison of the lattice 
quark mass function at various quark masses with 
a form derived from the Dyson-Schwinger equa- 
tion. The DSE describes the lattice data well 
at finite quark mass, but there is some differ- 
ence in the chiral limit. In particular the DSE 
solution changes more rapidly with momentum 
around one GeV. This suggests some non-linear 
behaviour at intermediate momenta and small 
quark masses, not captured by out simple extrap- 
olation. While such behaviour is certainly not 
surprising, it is also not well understood. 

4. ASYMPTOTIC BEHAVIOUR 

As we have already remarked, although the 
quark propagator is gauge dependent, and hence 
not observable, certain gauge invariant quantities 
can be extracted from it directly. In the contin- 
uum, the quark mass function has the asymptotic 
form [U, 



M(q 2 ) 



*oo q- 



[ln(g 2 /A 2 D )] 



(9) 



where fh is the RGI (renormalisation group in- 
variant) mass and the anomalous dimension of the 
quark mass is <1m — 12/(33 — 2Nf) for Nf quark 
flavours (zero here). In this form the fact that 
the mass function is independent of the choice of 
renormalisation point is manifest. The first term 
in Eq. is of nonperturbative origin and, in the 
chiral limit, a non-zero value for c indicates the 
spontaneous breaking of chiral symmetry. In this 
case it is related to the usual chiral condensate 
by 



[ln(/i 2 /A 2 



QCD 



)]' 



(10) 



where /z 2 is a choice of renormalisation point. The 
second term in Eq. |j5Jl comes from perturbation 
theory (one-loop) and breaks chiral symmetry ex- 
plicitly. The RGI quark mass in Eq. 10 can be 
replaced by the running quark mass, 



(11) 
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Table 2 

Results for the fit to Eq. ©. The fit region cor- 
responds to 1.9-2.9 GeV, which includes 51 data 
points, and was chosen to minimise the x 2 /dof. In 
the chiral limit this gives a value for the conden- 
sate of (-(i/n/')) 173 = 274(24) MeV at the reno- 
malisation point /i = 2 GeV. 



am 


a 3 c 


m (MeV) 


X 1 1 d.o.f. 





0.018(6) 





0.52 


0.012 


0.019(3) 


28(2) 


0.57 


0.018 


0.014(12) 


44(5) 


0.54 


0.024 


0.0096(47) 


59(3) 


0.52 


0.036 


0.0026(104) 


90(5) 


0.49 


0.048 


0.086(3) 


137(4) 


0.46 


0.072 


0.130(3) 


204(5) 


0.40 


0.108 


0.190(3) 


300(5) 


0.31 


0.144 


0.251(5) 


397(8) 


0.25 



Table [3 summarizes results for a fit of Eq. © to 
our lattice results. 

Having fit the asymptotic form to the quark 
mass function we should also be able to calculate 
the quark mass, but there is one piece of informa- 
tion missing from Tabled The effect of the bare 
quark mass on physical quantities is theory de- 
pendent; we have no a priori knowledge of where 
the physical point lies. We can match a physical 
quantity to its real world value and thus extract 
the relevant bare mass for that particular action, 
with those particular lattice parameters. In this 
case we choose the Asqtad quark action on the 
16 3 x 32 lattice at /3 = 4.60 and set the "physical" 
bare quark mass by extrapolation to the physical 
pion mass. 

We fit the five lowest ir masses to the form 

am?,- = V a 2 Bm 1 (12) 

the lowest order result from chiral perturbation 
theory. For the higher masses, the data devi- 
ates from this form due to the contribution from 
higher order terms; in the heavy quark limit cx 
to. The U(l) chiral symmetry of staggered quarks 
ensures that there is no additive mass renormali- 
sation. To find the physical point we find the bare 
mass corresponding to = 140 MeV for each 
jackknifc bin and get am = 0.00215(1). Then we 



fit the lowest four RGI masses in Table |21 to the 
form 

am = Aam (13) 

and find A = 1.57(16). Putting these together we 
get 

m = 5.3(6)MeV. 

The RGI quark mass in Eq. © is related to the 
running quark mass through Eq. i|ll|) . Thus by 
choosing the renormalisation point /i = 2 GeV 
and inserting A^ D = 239 MeV, we get 

m(fi = 2 GeV)^ = 3.1(4)MeV, 

or with Aq^ 1 = 691 MeV, 

m(p = 2 GeV) M ° M = 4.0(5)MeV. 

This result compares favourably with the Particle 
Data Group's number for the average of the up 
and down quark masses |19j . 

m(n = 2 GeV)^ = 2.5 - 5.5MeV. 

The curious behaviour of the condensate term 
at small masses is easily explained. The asymp- 
totic behaviour of the quark mass function 
quickly becomes dominated by the explicit chiral 
symmetry breaking mass; at the same time, large 
fluctuations make the mass function relatively 
noisy at these masses. Together, that means that 
sensitivity of the fit to the first term in Eq. (JJJJ is 
lost. At larger masses, the signal becomes cleaner 
and statistical errors shrink, restoring stability to 
the fits. One way of countering this effect would 
be to fix ac for each light mass to some reason- 
able value given the chiral and heavy quark val- 
ues. The difference to the end result would be 
small, but a future study should consider such 
systematics. 

This calculation is, of course, quenched. There 
are also systematic errors associated with the two 
extrapolations. Using other observables to set the 
scale - such as the p and nucleon masses, possibly 
supported by the appropriate finite-range regu- 
lated |20I21| staggered chiral perturbation theory 
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22 23] - would help us understand the system- 
atic uncertainty associated with the extrapola- 
tion. Such an approach could also provide a han- 
dle on the uncertainty associated with quenching. 
Finally, we should realise that this is a rather au- 
dacious approach: we are exploiting ultraviolet 
physics on a relatively coarse lattice. It is very en- 
couraging, however, that we are able to get such 
a promising result with such modest resources. 
The choices of improved quark and gauge actions 
are essential to this. 

5. CONCLUSIONS 

A comprehensive understanding of the 
quenched quark propagator in Landau gauge 
is now in place. Overlap and Improved Stag- 
gered actions produce propagators that are well 
behaved over a wide range of momenta. Despite 
this, more points in the deep infrared are re- 
quired to pin down precisely the behaviour near 
zero four-momentum. Also, the role of finite 
volume effects in the Z function at small quark 
masses needs further investigation. Calculations 
using dynamical configurations are underway. 
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